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We study theoretically the variation of the magnetic and the structural transition temperatures 
of the iron arsenide superconductors with external uniaxial stress. We show that the increase of the 
transition temperatures reported in recent experiments is compatible with a simple magneto-elastic 
model in which the physical origin of this variation is linked to the fact that antiferromagnetic bonds 
are longer than the ferromagnetic ones in the magnetic phase. We make predictions which can be 
verified in order to test the relevance of this model. 
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Introduction. At low doping, iron arsenide super- 
conductors exhibit magneto-structural transitions from 
paramagnetic metals with tetragonal crystalline symme- 
try to low temperature antiferromagnetic metals with 
(tt, 0) magnetic order and orthorhombic crystal structure 
(using the notation of IFe/cell Brillouin zone) [1]. These 
transitions are either concomitant or close to one another, 
thereby suggesting the presence of non-trivial magneto- 
elastic couphng [2, 3]. The structural transition consists 
of a C4-symmetry breaking that, from a purely electronic 
point of view, can be seen as a nematic transition [4]. 
Consistent with this interpretation, various experiments 
have reported in-plane anisotropy of the electronic prop- 
erties between the x- and the y- axes in the orthorhombic 
phase [5] . Since the study of such anisotropy is best per- 
formed using detwinned samples, experimental methods 
have been developed to perform the detwinning mechan- 
ically using in-plane uniaxial pressure [6]. However, an 
important prerequisite for interpreting these data is the 
understanding of how the magneto-structural transitions 
are themselves affected by the uniaxial pressure. Very 
recently experiments on BaFe2As2, one of the most ex- 
tensively studied systems for electron anisotropy, have 
shown that the both the structural and the magnetic 
transition temperatures, Ts and T/v respectively, are re- 
markably sensitive to uniaxial pressure [7, 8]. In fact, it 
has been reported that a modest compressive stress of 
0.7 MPa along the shorter y-axis results in an increase 
of Ts by about 10 K, and that of by few Kelvins [7]. 
Strikingly enough, tensile stress produces a similar in- 
crease of the magneto-structural transition temperatures 
[8]. This rules out a trivial magnetostriction effect be- 
cause, in that case, compressive and tensile stresses are 
expected to give opposite (increase vs. decrease) results. 

The purpose of this paper is to point out that these 
changes in Ts and Tpf can be understood quite readily 
as a consequence of magneto-elastic coupling. Specifi- 
cally, of the coupling that in the magnetic phase ensures 
that the antiferromagnetic bonds are longer compared 
to their ferromagnetic counterparts, a feature which is 



universal to all known iron based superconductors [1]. 
The magnitudes of the above changes are expected to in- 
crease if the system is near a second order orthorhom- 
bic instability, which is the case of BaFe2As2. As a 
corollary to these reasonings, we predict that applying 
tensile/compressive stress along the shorter/longer Fe-Fe 
bonds in the orthorhombic phase will, in fact, reduce the 
magneto-structural transition temperatures. 

Theory. We write the free energy of the system as 

F = Fm + Fe + Fme — ^xxUxx — ^yyUyyi 

where an denotes external uniaxial stress, with the con- 
vention that (7 < for compressive stress and cr > for 
tensile stress, and the components of the in-plane 

strain tensor with (i,j) — {x,y). The magnetic part of 
the free energy is given by 

Fm = A{\.\+\?_) I2 + --- , 

where L+ and L_ are magnetic order parameters cor- 
responding to (tt, 0) and (0, tt) orders respectively, and 
A = a{T — T^) where is the nominal Neel transition 
temperature. The elastic part is described by 

— 2 \^^^ ' ^yy } ' 2 v^^^ '^yy) 

~^ ^ \^xx ^^yy) ^ ' * * ; 

where K is the bulk modulus, Co is the temperature (T) 
dependent orthorhombic elastic constant, and i3 is a tem- 
perature independent constant. The important magneto- 
elastic part is given by 

Fme = ~g {uxx - Uyy) (L^ - L^) , (1) 

with the coupling constant g > such that in the mag- 
netic phase the antiferromagnetic bonds are longer than 
the ferromagnetic bonds (see below). The ellipses in the 
above equations denote terms that are not relevant for 
the current discussion. Since the systems of interest are 
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near a second-order structural transition where Co — > 0, 
magnetostriction effects are sub-dominant and can be ig- 
nored. 

The coupUng (1) has been shown to play a key role in 
establishing a universal phase diagram of these systems 
[2, 3]. In particular, the presence of tricritical points 
at which the magnetic transition changes character from 
first to second order were predicted for this phase dia- 
gram in [2], and it has later been confirmed experimen- 
tally in [9]. However, in the following we consider the 
case where the intrinsic instabilities (defined as the tem- 
peratures where the coefhcients A and Co are zero) are 
sufhciently apart in T such that the magnetic transition 
remains second order and the elastic response to exter- 
nal stress is linear at T/v- For the sake of concreteness 
we treat the case where > T/y (which is relevant for 
BaFe2As2). In the absence of external stress the mag- 
netic ordering takes place at 



Tjv(a = 0) = Tjv(0) 



n 



2gS/a, 



where S = {\Co\/By^-^ is the spontaneous orthorhom- 
bic distortion that appears for T < Tg due to the lat- 
tice instability (Cq < 0). Note that g > ensures that 
in the magnetic phase the antiferronagnetic bonds are 
longer than the ferromagnetic ones (mathematically, the 
> Uyy and the Uxx < Uyy solutions are coupled to the 



and L_ respectively). In the foUow- 



'^xx ^ ^yy 
instabilities of 
ing we describe two different experimental paths, namely 
stress applied above and below Tg (in the presence of 
finite stress T5 defines a cross-over where Co vanishes), 
and we assume that below Tg the system is in a single 
domain state with u^x > Uyy. 

(i) Uniaxial stress applied in the paramagnetic tetrago- 
nal phase. We consider uniaxial stress cr along one of the 
two equivalent directions, say y- for concreteness as in [7]. 
The strains induced by a are obtained by minimizing the 



total free energy 



OF 
dui 



0. Since we are assuming that 



the system first feels the underlying structural instability 
and then undergoes the magnetic transition (Tg > T/v), 
it is convenient to carry out this minimization by express- 



ing the total strain as Ui 



where 



repre- 



sent the strains that would appear spontaneously due to 
the structural transition at T5 (m°^ — Uyy = ±6). In 
this way, the change in the effective orthorhombic elas- 
tic constant Cq — > Cq = 2|Co| that takes place below 
Ts is captured via w?-. Sufficiently below T5 the elas- 
tic response is linear and the strains induced by a read 
Uyy = -{K + Co)/(if - Co)uxx = t{K + Co)/(4KCo). 
Substituting these expressions back in the free energy 
we find that the effective magnetic energy is modified to 
F'j^ = Fm + ffcr(L^ - L?.)/(2Co). The Neel temperature 
is determined by the leading instability between L+ and 
L„ associated to (tt, 0) and (0, tt) orders respectively. In 
the case of compressive stress the (tt, O) magnetic state is 
promoted, while with tensile stress it is the (0, tt). But for 
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uniaxial stress (along short Fe-Fe bonds) 

FIG. 1: (Colour online). Schematic variation of mag- 
netic and structural transition temperatures, Tn (red) and 
Ts (blue) respectively, with uniaxial stress in detwinned sam- 
ples. According to Eq. (2) Tn increases, and Ts tracks it 
(solid lines), for both compressive and tensile stress applied 
in the tetragonal phase as seen experimentally [7, 8] . Eq. (3) 
predicts that Tn decreases, and Ts tracks it (dash lines), for 
tensile/compressive stress appHed along the short /long Fe-Fe 
bonds in the orthorhombic phase. 



both compressive as well as for tensile stress increases 
as 



Tn{<j)^Tm{0) 



Co a 



(2) 



Note that this conclusion is also valid for systems where 
the intrinsic Neel transition appears before the structural 
one (i.e., T% > Ts), with Co replaced by Co in Eq. (2). 

(ii) Uniaxial stress applied in the paramagnetic or- 
thorhombic phase. In the single-domain state with Uxx > 
Uyy, which promotes the (tt, 0) magnetic order, the uniax- 
ial stress can be applied either along the long or the short 
directions {x and y respectively). The induced strains 
can be obtained by minimizing the total free energy as 
before. After substuting the resulting expressions back 
in the free energy we find that stress applied along the 
short y direction shifts the Neel temperature to 



TN{a)=TNiO) 



Co a 



(3) 



Thus, while T^ increases in the case of compressive 
stress, it actually decreases in the case of tensile stress 
unlike in the previous case (i). Obviously, the conclu- 
sion is reversed if the stress is applied along the long 
x-direction. 

Discussion. Eqs. (2) and (3) are the main results of this 
paper, which are illustrated in Fig. 1. Note that these 
results are a direct consequence of the magnetoelastic 
coupling described by (1) and symmetry considerations. 
As such, they do not depend upon whether the magnetic 
excitations are more itinerant-like or more localized-like. 
Furthermore, it can be argued that these conclusions re- 
main unchanged even if the C4-symmetry breaking is due 
to electron spin nematicity or due to orbital ordering (in 
which case orthorhombicity is a secondary order param- 
eter). The linear increase of T/v with both compressive 
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uniaxial stress 



FIG. 2; (Colour Online). Schematic variation of the mag- 
netic transition temperature Tjv with uniaxial stress a for 
different dopings. The slope dTN/da, which is inversely pro- 
portional to \Ts — Tm\ in the current theory, is predicted to 
decrease with doping. 



and tensile external stress described by Eq. (2) has been 
observed experimentally for small a [7, 8] . Note that the 
sign of g is crucial for understanding these experimental 
findings. The decrease of Tjy with compressive/tensile 
uniaxial strain applied along the long/short Fe-Fe bonds 
described by (3), on the contrary, is a prediction to be 
confirmed in future experiments. In fact, in the case 
where the system is not completely detwinned, the ap- 
plied stress reduces the Neel temperature of the "wrong" 
domains compared to that of the "correct" ones, and this 
results in the smearing of the experimentally determined 
Tjv. Furthermore, since the orthorhombic elastic con- 
stant Co appears in the denominator in both Eqs. (2) 
and (3), closer the system is to a second order orthorhom- 
bic instability as in BaFe2As2 and EuFe2As2 [10], larger 
is the variation of T/y with a. On the contrary, we ex- 
pect the magneto-structural transitions to be relatively 
insensitive to a for CaFe2As2 and SrFe2As2 for which the 
structural instability is first order. In fact, if we assume 
that at the Neel transition Cq varies linearly with tem- 
perature as Co oc T — Ts, then we get that the slope 
dT^ / da varies inversely with (T5 — Tjy). This conclusion 
can in principle be tested by comparing the results for 
samples at various doping as illustrated in Fig. 2. 

Finally, we focus on the variation of Tg with exter- 
nal uniaxial stress. It has been pointed out that the 
orthorhombic transition can be interpreted as a result 
of structural softening due to magnetic fluctuations that 
develop near a (tt, 0) magnetic order [11, 12]. This in- 
terpretation provides a simple explanation why T5 and 
T/v track each other in the temperature-doping phase 
diagram of all the iron arsenides. More concretely, if 
we assume that the structural transition occurs in a 
regime where the magnetic fluctuations are controlled 
by the magnetic quantum critical point, we get that the 
correction to the orthorhombic elastic constant is [11] 
5Co = -AIii[Tq/{T - Tn)], where A and Tq are con- 
stants that depend on microscopic details. The condition 



Co + SCq = determines the structural transition point, 
from which we obtain 

Ts =Tw+roe-^«/^. (4) 

Thus, we expect Ts to increase (or decrease) with a in 
order to keep track of the increase (or decrease) of T^r. 
Note that this conclusion is qualitatively valid even if the 
magnetic fluctuations are controlled by the finite temper- 
ature magnetic critical point. 

Conclusions. We have studied the variation of the 
magneto-structural transition temperatures of the iron 
arsenide systems with uniaxial stress. We have shown 
that the physical origin of this variation is the magneto- 
elastic coupling which makes the antiferromagnetic bonds 
longer than the ferromagnetic ones in the magnetic 
phase. The conclusions of the theory are compatible 
with the current experimental findings. The magnitude 
of this variation is relatively large because of the prox- 
imity of the orthorhombic instabilty, which is the case of 
BaFe2As2 (and presumably also of EuFe2As2). In or- 
der to test this interpretation we predict that (i) the 
transition temperatures will decrease if, after preparing 
the sample in a single-domain state, tensile/compressive 
stress is applied along the shorter/longer Fe-Fe bonds in 
the orthorhombic phase, (ii) and that the slope dT^/da 
is inversely proportional to \Ts — Tjvl when the two tran- 
sitions are close enough (but not too close to make the 
magnetic transition first order) . Note that our results are 
independent of microscopic details concerning the mag- 
netic properties of the iron arsenides. 

We thank Y. Gallais, G. Garbarino and E. Kats for 
fruitful discussions. 

Note added. Hu et al. [13] have also studied pressure 
effects in iron-pnictides from a microscopic spin model 
and a pure electron-nematic point of view. Their results 
are compatible with Eq. (2). 
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